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ABSTRACT 



Context. Determination of ages of powerful radio sources is crucial in understanding galaxy evolution, the activity cycle of 
galactic nuclei and their influence on the surrounding intergalactic medium. So far, several different methods for age estimates 
of classical double radio galaxies have been proposed and widely used in the literature, although each of them faces some 
difficulties due to observational limitations and/or freedom in choosing the underlying model assumptions. 
' Aims. We propose a new approach in determining ages of FR II type radio sources, on one hand exploiting a dynamical model 
f*") , developed for these objects by Kaiser et al. (1997), and on the other hand using multifrequency radio observations not necessarily 
restricted to the high-resolution ones. 

Methods. In particular, we apply the assumed dynamical model to a number of FR II type radio galaxies observed at different 
radio frequencies, and fit — for each frequency separately — the model free parameters to the observed sources' quantities. Such 
Q ' a procedure, using enlarged in fact a number of observables, enables us to determine relatively precise ages and other crucial 
S—j \ characteristics (like the jets' kinetic power) of the analyzed sources. 

(/} . Results. The resulting age estimates agree very well with those obtained by means of 'classical' spectral ageing method for objects 
' not older than 10 Myr, for which good-quality spectral data are available. The presented method is however also applicable in 
the case of the sources older than this, and/or the ones for which the only available low-resolution radio data do not allow for 
? T ~j ■ detailed spectral ageing studies. Interestingly, the estimated ages always correspond to the realistic values of the jets' advance 
[ velocity ~ 0.01 c — 0.1 c 

Conclusions. Our analysis indicates that the main factor precluding precise age determination for FR II type radio galaxies 
regards the poorly known shape of the initial electron energy distribution injected by the jet terminal shocks to the expanding 
lobes/cocoons. We briefly consider this issue, and conclude that the broad-band single power-law form assumed here may be 
accurate enough for the presented estimates, although most likely it does not strictly correspond to some well-defined realistic 
particle acceleration process. Instead, it should be considered as a simplest model approximation of the initial electron continuum, 
averaged over a very broad energy range and over the age of the source, with the effective spectral index which may be different 
for different sources, however within the relatively narrow range p = 2.0 — 2.4 suggested by our modeling. 

Key words, galaxies: active - galaxies: evolution - radio continuum: galaxies 



1. Introduction 

The dynamical age of extragalactic radio sources is of 
a great importance for studying and understanding de- 
tails of the physical processes taking place in both Active 
Galactic Nuclei (AGNs) emitting energy in forms of jets 
boring through the ambient intergalactic medium (IGM), 
and also lobes or cocoon consisting of shocked (backfiow- 
ing) jet material and shocked IGM. The 'standard model' 
of classical double radio sources (Blandford & Rees 1974; 
Scheuer 1974) implies that the axial lengthening of the co- 
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coon is governed by the balance between the jet's thrust 
and the ram pressure of the IGM, whilst the width of the 
cocoon is determined by its internal pressure. This sce- 
nario has formed the basis for all further analytical models 
of the dynamics and radio-emission properties of powerful 
double radio sources (e.g., Begelman & Cioffi 1989; Falle 
1991; Nath 1995; Kaiser et al. 1997; Chyzy 1997; Blundell 
et al. 1999; Manolakou & Kirk 2002, Kino & Kawakatu 
2005). 

All such dynamical models allow one to predict a num- 
ber of observational parameters of the sources (e.g., linear 
size, width of the cocoon, radio power at different wave- 
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lengths) for their given ages and redshifts, assuming val- 
ues of the models' basic free parameters (e.g., jet power, 
external medium density and its gradient, etc.). A com- 
parison of the predictions of the three main scenarios, i.e. 
the ones by Kaiser et al. (1997), Blundell et al. (1999) 
and Manolakou & Kirk (2002), was recently undertaken 
by Barai & Wiita (2005). They argued that no existing 
model could give acceptable fits to the common proper- 
ties of radio sources from the survey considered (3CRR, 
6CE, 7CRS), though the Manolakou & Kirk model gave 
'better overall results than do either of the other two'. 

The age of a classical double radio source can, in prin- 
ciple, be determined from its apparent radio spectrum. 
There is no doubt that a radio continuum spectrum in 
different parts of a radio source should contain crucial in- 
formation about its history, depending on various energy 
losses (and gains) of the radiating particles. This simple 
idea, relying on the determination of a single characteris- 
tic break frequency in the radio spectrum (at a given part 
of the source) caused by the above time-dependent losses, 
was used to estimate the radiative ages and expansion 
speeds in large samples of powerful 3CR sources (Myers 
& Spangler 1985; Alexander & Leahy 1987; Leahy et al. 
1989; Carilli et al. 1991; Liu et al. 1992), as well as in sam- 
ples of low-luminosity and intermediate-luminosity radio 
galaxies (Klein et al. 1995; Parma et al. 1999). However, 
the spectral break and steepening of the spectrum beyond 
this break need not be entirely due to radiative ageing. 
A possible role of the local magnetic field structure and 
evolution, details of the backflow process of the cocoon 
material, or the difficulties in disentangling various en- 
ergy losses of the radiating particles have been raised in a 
number of papers (e.g. Rudnick et al. 1994; Eilek & Arendt 
1996; Jones et al. 1999). 

The task of resolving some of the above difficulties 
was undertaken by Kaiser (2000; hereafter referred to as 
K2000) who has constructed a 3-dimensional model of the 
synchrotron emissivity for a cocoon of a FR II-type radio 
source, and applied it to the self-similar model of the co- 
coon's dynamics developed by Kaiser & Alexander (1997; 
hereafter referred to as KA) and its extension by Kaiser 
et at. (1997; hereafter KDA). In the K2000 model, the 
cocoon is split into small volume elements whose surface- 
brightness evolution is traced individually. The author has 
argued that the proposed careful analysis of the bulk back- 
flow process and of the electron energy losses (both radia- 
tive and adiabatic) can provide an accurate estimate of 
the sources' spectral ages, being in addition in agreement 
with the estimates of their dynamical ages. Blundell & 
Rawlings (2000) have contended however that this could 
only be the case if these ages were much less than 10 Myr. 

Analyzing the K2000 model, we realized that in the 
majority of the extended FR II-type radio sources — even 
those without distorted lobe structures — the surface- 
brightness profiles are far from the expected smooth 
shapes, causing the fitted free parameters to be highly 
uncertain. In addition, the K2000 method requires rather 
high resolution observations of the radio lobes. But usu- 



ally high-resolution observations of low-brightness sources 
cause a serious loss of the flux density, so that application 
of the K2000 method is in fact limited to a few strongest 
sources such as Cyg A. 

The KDA model has been successively used by 
Machalski et al. (2004a, b) to derive the basic physical 
parameters of the FR II-type radio sources with known 
age, namely the jet power, the central density of the am- 
bient medium, the cocoon's energy density and pressure, 
and the total energy deposited in the radio source. All of 
these values were derived by means of fitting the model 
free parameters to the observed parameters of the source: 
its redshift, monochromatic radio luminosity at 1.4 GHz, 
projected linear size, and the cocoon's axial ratio. In this 
paper we propose a simplified method of the age determi- 
nation for powerful radio sources based on a further ex- 
ploitation of the KDA model, which avoids the limitations 
present in the K2000 method. In particular, we show that 
performing the above fitting procedure for a few observing 
frequencies, i.e. enlarging the number of the observables 
used by adding the radio spectrum data of the entire co- 
coon, a unique set of the age and the values for the two 
fundamental model parameters — the jet power and the 
central density of the ambient medium — can be found. 
This, in turn, enables us to determine relatively precisely 
the age of the analyzed source, even if the available (low- 
resolution) radio data do not allow for detailed spectral 
ageing studies. 

The paper is organized as follows. First, a brief pre- 
sentation of the analytical KDA model together with its 
application to the dynamical ageing studies are given in 
Section 2. In Section 3 we present the dataset used in the 
subsequent analysis. Exemplary results of the dynamical 
age determination for the radio galaxy Cygnus A and its 
dependence on the assumed model parameters are given 
in Section 4. Application of the proposed method to a va- 
riety of other FR II type radio sources is summarized in 
Section 5. Remarks on the implications of the obtained re- 
sults for modeling the electron injection spectrum in lobes 
of powerful radio sources are given in Section 6. Final con- 
clusions are presented in Section 7. 



2. The KDA model and its application 

2.1. Source dynamics 

The dynamical expansion of a FR II-type radio source 
built in the KDA model is based on the dynamics de- 
scribed in the self-similar KA model. Based on the original 
idea of Blandford & Rees (1974) and Scheuer (1974), it as- 
sumes that the classical double radio structure is formed 
by twin jets emerging from the AGN in opposite directions 
into a surrounding tenuous environment. The jets termi- 
nate in strong shocks where the jet particles are shocked, 
accelerated, and finally inflate the cocoon. A density dis- 
tribution of the unperturbed external gas is approximated 
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by a power-law radial distribution scaling with distance 
r > a from the center of the host galaxy as 



P(r) = Po 



0() 



where po is the central density at the core radius ao, and 
the exponent describes the density profile in the sim- 
plified King's (1972) model. Such a distribution of the 
ambient medium is assumed to be invariant with redshift. 

Half of the cocoon (from the core to the hotspot) is 
approximated by a cylinder of length Lj et — D s /2 and 
base diameter b. Here D s is the total deprojected length 
(linear size) of the source. The ratio of the source's length 
to the base diameter is called the 'axial ratio' AR; in the 
original KDA paper, half of the cocoon has axial ratio of 
i?T = AR/2. The cocoon expands along the jet axis driven 
by the hotspot plasma pressure ph and in the perpendicu- 
lar direction by the cocoon pressure p c . The rate at which 
energy is transported from the AGN along each jet, i.e. 
the jet power Qj et , is assumed to be constant during the 
source's lifetime. The model predicts self-similar expan- 
sion of the cocoon and gives analytical formulae for the 
time evolution of various geometrical and physical param- 
eters, e.g. the length of the jet (Eqs. 4 and 5 in KA) 



1/(5-/3) 



£3/(5-/3) ^ 



(1) 



and the cocoon pressure (Eq. 20 in KA) 



Pc = 



18 c 



(2-/3) 



(r x + i)(5-/?) 2 7? h( 



(poao) 



3/(5-/3) 



(4+/3)/(5-/3) 



(2) 



where ci is a dimensionless constant (Eq. 25 in KA), T x is 
the adiabatic index of the unshocked medium surrounding 
the cocoon, t is the time elapsed since the jet had started 
from the AGN, i.e. it is the source actual age, and The = 
Ph/Pc is the pressure ratio. In our calculations we use the 
empirical formula taken from K2000 

P hc = (2.14 - 0.52/3)^-04-0.25/3^ 



2.2. Source energetics and radio power 

The energy delivered by the twin jets during the source 
lifetime is simply equal to 2Q iet t. Accounting for the en- 
ergy losses of relativistic particles due to adiabatic ex- 
pansion within the time-dependent volume of the cocoon 
V c (t), the ratio of the jets' energy and the energy stored 
in the cocoon E to t is 



2Qj Ct t _ 9T C - 4 - 9(T C - 1) 
E tot 5-0 4(5-/3) 



Phc, 



where E tot (t) = u c (t)V c (t), u c = p c {t)/{T c - 1) is the 
energy density in the cocoon, and T c is the adiabatic index 
of the cocoon's matter. 



Following KDA, we assume that the cocoon is com- 
posed by many small volume elements filled with mag- 
netized plasma and cosmic rays which are injected into 
the cocoon at the jet head — the jet terminal shock. The 
initial energy distribution n(7i) of injected particles with 
the Lorentz factor 7; is assumed to be of a single power- 



law form with the spectral index p: n(7i) 



"o7i 



The 



unknown normalization parameter no is given by integrat- 
ing the initial power-law distribution over the entire (as- 
sumed) energy range [7 

,min, 7i,maxJ ■ 

The relativistic particles injected into the cocoon loose 
their energies in a magnetic field of energy density ub, 
and thus the cocoon's radio spectrum evolves with time. 
Tracing the effects of adiabatic expansion, synchrotron 
losses (with the assumed effective isotropisation of the par- 
ticles' pitch angle distribution), and inverse Compton scat- 
tering on the cosmic background radiation in the volume 
elements independently, the radio power of the cocoon P v 
at a fixed observing frequency v is obtained by summing 
up the contributions from all volume elements, resulting 
in an integral over time (see Eq. 16 in KDA). It depends 
on the source's age t and redshift z, the jet power Qj et , 
the half-cocoon axial ratio i?T, the exponent in the initial 
power-law distribution of relativistic particles p, the val- 
ues of the low-energy and high-energy cut-offs 7i jm i n and 
7i.max, and on the ratio £ of the magnetic field energy ub 
to the energy of relativistic electrons u e . The integral is 
not analytically solvable and has to be calculated numer- 
ically. 

On the basis of the above analytical model we aim to 
predict the dynamical age and the two physical param- 
eters of the analyzed FR II-type radio sources: Qj et and 
Po- In order to derive these parameters, all the other free 
parameters of the model (r c , T x , a , 0, C, P, 7i,min, 7i,max) 
have to be fixed (see Table 1). 

2.3. Selection of the free model parameters 

Following the KDA analysis, we adopt their Case 3 where 
both the cocoon and the ambient medium are assumed to 
be best described by the 'cold' equation of state (r c = 
T x = 5/3). For the initial ratio of the energy densities 
of the cocoon's magnetic field and particles we use the 
equipartition condition ( = ub/u c = (1 +p)/4, well sup- 
ported by the X-ray observations of the lobes in powerful 
radio sources (see Kataoka & Stawarz 2005, Croston et al. 
2005, and references therein). 

Selection of the particular (universal) values regard- 
ing external medium parameters ao and requires some 
justification. We note, that even careful 2-D modeling of 
a distribution of radio emission for powerful sources with 
quite regular structure can lead to values of ao discrepant 
with those implied by the X-ray observations of intergalac- 
tic thermal environment. In our simplified approach we 
assume ao = 10 kpc for all the analyzed sources, a con- 
servative value between 2 kpc used by KDA and 50 kpc 
found by Wellman ct al. (1997). Wc also use a constant 
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value of p = 1.5, in agreement with Daly (1995), Blundcll 
et al. (1999), Willott et al. (1999), and K2000. Note, that 
the assumed ambient medium profile is then slightly flat- 
ter than the one adopted in the original KDA paper on 
the basis of Canizares et al.'s (1987) analysis, who found 
ft = 1.9 to be typical at about 100 kpc from the galactic 
nucleus (but see Gopal-Krishna & Wiita 1987). 

Another free parameter of the model is the orientation 
of the jet axis to the observer's line of sight 9. We assume 
9 = 90° for the giant-size radio galaxies and 9 = 70° for 
'normal'-size radio galaxies. This latter value is justified 
on the basis of the unified scheme for extragalactic radio 
sources. In this scheme, the average orientation angle for 
radio galaxies is (0-rg) — 69° (Barthel 1989). The appar- 
ent size D of a radio source then yields the model cocoon 
length 

Ljet =£>/(2sin0). (3) 

The next parameters which significantly influence the 
model predictions (as will be shown and discussed in the 
next sections) are the initial power-law exponent p, as well 
as the low- and high-energy cut-offs of the initial electron 
energy distribution, y hmin and 7;, max . 

The cut-off energies in the electron energy distribution 
injected by the jet terminal shock to the expanding cocoon 
are poorly constrained by both observations and theory. 
The values considered in the literature for 7L max vary from 
10 4 up to > 3 x 10 8 (see Barai & Wiita 2005). Here we 
initially study a very broad range 7i imax = 10 4 — 10 10 (sec- 
tion 4), in order to check how the results of the proposed 
method depend on the choice of this particular model pa- 
rameter. We found that the effect is minor. We note in this 
context that the high values of the maximum Lorentz fac- 
tors for electrons accelerated at the jet terminal shocks are 
more likely, since a number of hotspots in powerful radio 
sources are established sources of synchrotron X-ray emis- 
sion, and hence since production of high-energy (ji > 10 7 ) 
electrons thereby seems to be a general property of radio 
galaxies and quasars (see discussion in Hardcastle et al. 
2004, Kataoka & Stawarz 2005, and references therein). 
For these reasons, later on in the analysis (section 5) we 
consider a slightly narrower range of 7;, max = 10 7 - io 8 - 5 . 

As for the minimum electron energy, we note that most 
authors assume 7i. m i n — > 1, although such a choice is 
made largely for simplicity reasons. Interestingly, recently 
Blundcll et al. (2006) argued that the X-ray observations 
of giant radio galaxy 4C 39.24 indicate an initial electron 
low-energy cut-off as high as 7;, m i n ~ 10 4 (see also in this 
context Carilli et al. 1991 for the low-frequency radio ob- 
servations of CygA hotspots). We believe, however, that 
these data do not show unambiguously any sharp low- 
energy cut-off of the electron energy distribution at some 
particular (ultrarclativistic) energy, but instead indicate 
only spectral flattening occurring at low (7i jm i n < 10 4 ) 
electron energies. Thus, we take below 7i, m ; n = 1, assum- 
ing that the particles at the jet terminal shock are picked 
up by the acceleration process(es) directly from the 'ther- 
mal pool' of the cold plasma carrying the bulk of the jet 



kinetic energy, to form the non-thermal electron energy 
distribution injected further to the expanding cocoon. We 
checked however that taking 7i jm i n = 10 (as favored by 
the analysis of Barai & Witta 2005) does not change our 
results obtained with 7i im i n = 1 presented in the next sec- 
tions. 

Table 1. Observational and model parameters (see sec- 
tions 2.1-2.3 for detailed descriptions). 





Symbol 


Dimension 


Observational 


2 


[dimensionless] 


parameters 


D 


[kpc] 


derived 


AR 


[dimensionless] 


from radio maps 


V 


[MHz] 


and spectrum 




[WHz-'sr- 1 ] 


Model free 


a 


[kpc] 


parameters 


P 


[dimensionless] 


(to be fixed) 


p=2a inj + l 


[dimensionless] 






[dimensionless] 






[dimensionless] 




c 


[dimensionless] 




t 


[yr] 




e 


[deg] 


Model parameters 


Qict 


[W] 


derived for given 


Pa 


[kgm 3 ] 


values of t, 


Pc 


[Nm- 2 ] 


ainj and 7i, ma x 


u c 


[Jm- 3 ] 




Etot 


[J] 



The above choices regarding low- and high-energy elec- 
tron cut-offs seem to be justified, although they may lead 
to inconsistency with the main assumption of the model, 
namely with the assumed single power-law form of the ini- 
tial electron energy distribution: if the low-energy spectral 
flattening of the electron spectrum is the case (as has to 
be accepted if indeed 7i jTO i„ = 1), it suggests that the 
particle distribution injected by the jet terminal shocks to 
the lobes/cocoons is characterized by some more complex 
spectral shape, for which a broken power-law could be for 
example a simple approximation. Meanwhile, in the origi- 
nal KDA model, universal single power-law initial energy 
distribution of the relativistic particles with the energy 
index p = 2.14 was assumed, corresponding to the ra- 
dio spectral index a- m j = 0.57. This value agrees perfectly 
with the one obtained for the classical FR II-type radio 
source Cyg A by means of fitting its radio spectrum by the 
'continuum injection' (C.I.) model (for the spectral age- 
ing theory review see Myers & Spangler 1985; Carilli et 
al. 1991). It should be noticed that p = 2.14 is close to, 
although not exactly equal to, the model spectral index of 
relativistic particles accelerated at non-relativistic shocks 
by the 1st order Fermi (diffusive) acceleration process. 

Most importantly, the low-frequency radio observa- 
tions of a number of powerful radio sources imply a large 
spread of the injection spectral index, 0.38 < a; n j < 0.95. 
This indicates that there may be no universal value of this 
parameter for all cocoons. Moreover, we note that the low- 
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frequency radio spectra usually steepen with redshift (see 
Laing & Peacock 1980; Macklin 1982; Gopal-Krishna 1988; 
van Brcugel & McCarthy 1990). The physical reason for 
such a dependence is unclear, although the effect seems to 
be real. Therefore, keeping in mind all the uncertainties 
regarding ai n j (or equivalently p) which reflect our poor 
understanding of the particle acceleration processes tak- 
ing place at mildly relativistic collisionless astrophysical 
shocks, and observational limitations, we decided to treat 
the injection spectral index as a free parameter, which 
may be different for different sources. 

We also emphasize, that the particular value of power- 
law index p obtained for a given source by the fitting 
procedure does not correspond necessarily to some well- 
defined acceleration process taking place at the jet termi- 
nal shock. Instead, it should be considered as an effective 
spectral index of the injected electron continuum (char- 
acterized most likely by more complex spectral shape), 
averaged over a very broad energy range and over the life- 
time of the source. This issue will be discussed further in 
Section 6. 

2.4. The fitting procedure 

In Table 1 we summarize all the observational and model 
parameters characterizing the source's cocoon and used in 
the present calculations. 

Assuming a value of the source age t, and having fixed 
all the free parameters of the model as discussed in section 
2.3, we find the jet power Qj t and the central density of 
the external medium po for an investigated source by an 
iterative solution of the system of two equations: (i) Eq. 1 
equated to Eq. 3 for the jet length Lj et , and (ii) the inte- 
gral for the luminosity of the cocoon P„ (Eq. 16 in KDA, 
see Sect. 2.2), requiring the matches of the solution to the 
observed values of D and P v at a number of different fre- 
quencies, respectively. The above fitting procedure proved 
to be stable and always provides unique solution. 

Performing the above calculations 1 for a number of 
age values t assumed for a given source, we obtain a set of 
solutions for Qj c t(£) and po(t), which show, as expected, 
a correlation of Qj e t and anticorrelation of po with the 
source age t for a given apparent size D and total co- 
coon luminosity P v . An example of this effect for the four 
observing frequencies is shown in Fig. 1. A perfect inter- 
section of all the Qj e t — po curves corresponding to dif- 
ferent observing frequencies at some particular 'true' age 
of a source is expected if the observed shape of its total 
radio spectrum (used here to determine the radio power 
at different frequencies) agrees with the theoretically pre- 
dicted one in the framework of the C.I. model, and if the 
dynamical (KDA) model applied here is indeed correct. 



1 Our programme DYNAGE to calculate the model param- 
eters (Qjet, Po, Pc, Uc, and -Etot) for a set of given values 
of t, Qi n j, and 7i im ax as described in this paper and illus- 
trated in Fig. 1, is available upon request from J. Machalski 
(machalsk@oa.uj .edu.pl). 



However, the observed spectrum of a number of the ana- 
lyzed sources evidently departs from the predicted behav- 
ior. In such cases, the sources' high-frequency spectra are 
in particular steeper than expected, and hence the resul- 
tant Qjot - Po curve corresponding to the highest observing 
frequency bisects other curves at the age inconsistent with 
that implied by the remaining ones. 



3C330 



11 



40 



Of 



o 



39 




151 MHz 
408 MHz 
♦ 1400 MHz 
□ 4800 MHz 



38 



-24 -22 -20 

log(p [kg/m 3 ]) 



-18 



Fig. 1. A set of Qj C t(t) and po(t) solutions at the four 
observing frequencies, for a number of the age values as- 
sumed for the radio galaxy 3C 330. 

A 'goodness' of the intersection is quantified by the A 
measure which is defined as follows: indicating logpo(i) = 
x(t) and logQjot(t) = y(t), we calculate the mean sep- 
aration between the points determined by the values of 
[x(t),y(t)] in the log<5j ct -logp plane as 



A(i) 



n(n — 1) 



n— 1 / n— 1 



:J2 £{M*) - *,+iW] 2 + [Viit) V 3+ i{t)] 2 } 



1/2 



,(4) 



where n is the number of the observing frequencies in- 
cluded. A minimum of A(t) (i.e. 'the best' intersection of 
the Qjet — Pa curves for different frequencies) is consid- 
ered as an indicator of the source's real dynamical age t. 
This minimum also discriminates a unique values of <5j t 
and po, as well as p c and i?tot for the considered source. 
Hereafter, the age corresponding to such a set of Qj et and 
po values is called the age solution. 

However, the calculations indicate that the 'age solu- 
tion' significantly depends on the exponent in the initial 
power-law energy distribution, p. In addition, they suggest 
that there is no canonical value of this injection spectral 
index which would be good for each and every source. 
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Table 2. List of sources and their observational parameters derived from radio maps and spectra. The radio power is 
in units of W Hz -1 sr _1 



Source name z D[kpc] AR g[°] logPisi jogjog logPi4uo jog-FW) 



CygA 


0.0564 


135 


3.9 


80 


27.796 


27.465 


26.963 


26.358 


3C55 


0.7348 


496 


6.4 


70 


27.644 


27.247 


26.718 


26.148 


3C103 


0.330 


458 


7.0 


70 


26.906 


26.630 


26.140 


25.552 


3C165 


0.296 


384 


3.4 


70 


26.579 


26.244 


25.769 


25.227 


3C239 


1.786 


96 


2.7 


70 


28.483 


28.092 


27.541 


26.914 


3C247 


0.749 


95 


4.8 


70 


27.350 


27.117 


26.746 


26.284 


3C280 


0.996 


113 


2.9 


70 


27.990 


27.697 


27.288 


26.829 


3C292 


0.71 


960 


8.4 


90 


27.406 


27.020 


26.547 


26.075 


3C294 


1.779 


135 


3.8 


70 


28.412 


27.984 


27.431 


26.841 


3C322 


1.681 


283 


5.0 


70 


28.171 


27.906 


27.485 


26.954 


3C330 


0.549 


395 


10.9 


90 


27.432 


27.150 


26.771 


26.334 


3C332 


0.1515 


229 


4.9 


70 


25.744 


25.460 
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(This in fact agrees with the preliminary notes given at 
the end of the previous subsection 2.3; this issue will be 
discussed further in section 6.) For this reason, the best 
value of p (or equivalently cti n j), needed for a precise esti- 
mate of the dynamical age of a particular FR II-type radio 
source, is determined below by finding the minimum of the 
jet kinetic energy delivered to the cocoon (see section 4.3). 
The age solution corresponding to the minimum jet energy 
is called hereafter the best solution. 



3. The data 

The data necessary for the age determination of a given 
radio source with the method presented above are: its pro- 
jected linear length D, its axial ratio AR, and its total 
radio power P v at a number of observing frequencies, and 
the source's redshift, z. 

As we limit our method to FR II-type sources only, D 
is invariant of frequency and can be precisely determined 
from the radio maps. The lateral size (base diameter b of 
the cocoon) is taken as the average of the widths of the two 
lobes. They are determined on the maps with the highest 
sensitivity to surface brightness as the deconvolved width 
of a transversal cross-section through the lobes measured 
half-way between the core and the outer edge of the source. 
Thus, AR = D/b. 

The source power P v is determined at four frequencies: 
z/=151, 408, 1400, and 4800 MHz. Because not all sources 
were observed at these particular frequencies, the cocoon's 
spectrum has been derived from the available flux densi- 
ties with the contributions from the radio core and the 
hotspots subtracted if possible. Next, the simple analyti- 
cal functions: y = a+bx+c exp(±x), where x = log ^[GHz] 
and y = log S(v) [mJy], have been fitted to these flux den- 
sities weighted by their given error. These two functions 
were found to provide the best fit (the least x 2 ) to the 
data. Then, the values of P151, -P408, fi400 and P4800 have 



been calculated using flux densities found with the above 
functions and the appropriate radio K-corrections 2 . 

The FR II-type radio sources analyzed in this pa- 
per, and their observational parameters as used in the 
fitting procedure, are given in Table 2. These sources, 
subjectively selected from the sample of Machalski, 
Chyzy & Jamrozy (2004a), include high- redshift and low- 
redshift, high- luminosity and low- luminosity, giant-sized 
and normal-sized radio galaxies. The sources in Table 2 
do not include quasars whose jets' axis (and hence radio 
structure) may be strongly projected onto the sky. 

4. Fitting result for CygA 

In this section the dynamical age estimate and its depen- 
dence on the values of the model free parameters is pre- 
sented for CygA. 

4.1. The age solution and its dependence on a- m j and 

Ti,max 

Assuming 9 = 80° and adopting (after K2000) parameters 
p = 2.14 and 7i !max = 10 4 , the fitted values of Qj et and 
po for different hypothetical ages of this source are shown 
as the log Qjet-log po diagram in Fig. 2. This figure shows 
that the four curves corresponding to the four observing 
frequencies intersect at the age of about 10 Myr. The age 
solution giving a minimum of the A(t) measure for the 
set of the model free parameters (recalled in the caption 
of Fig. 2) is t = 10.9 ± 0.9 Myr. 

As mentioned previously, a minimum of A(i) appears 
to be dependent on the assumed parameters of the ini- 
tial energy distribution of relativistic particles. The de- 
pendence of A on t for two different (rather extreme) val- 
ues of 7i, m ax = 10 4 and = 10 10 is plotted in Fig. 3 for 

2 There is no need to use the same frequencies for different 
sources; a few observing frequencies should cover a wide range 
of the radio spectrum 
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Fig. 2. Qj ot — p diagram for CygA at the four observ- 
ing frequencies. The dotted diagonal lines indicate an age 
related to the fitted Qj C t and po values. Other free param- 
eters for all the curves are: Rt—1.95, (3=1.5, aii n j=0.57, 
and 7; max =10 4 (see also section 2.3) 

different values of a; n j . As shown, the age solution is not 
very sensitive for the particular choice of 7i jma x, especially 
if Ctinj > 0.5. 

The dependence of the resulting age solution on a; n j for 
the same two extreme values of 7i jmax is shown in Fig. 4. 
It is clear that the age solution strongly increases with 
the decreasing a; n j , and hence that the minimum A mea- 
sure is smaller for larger assumed values of ai n j. This is 
due to the fact that in the case of flatter and flatter ini- 
tial electron energy distributions more and more time has 
to elapse for the radio continuum to steepen up to some 
given (observed) level. As a result, the obtained age solu- 
tion anticorrelates with the injection spectral index, and 
hence also the minimum of A decreases with decreasing t. 
Such a decrease continues until reaching the point where 
the radio continuum did not have time to steepen at all, 
correspondingly to the situation when all the curves on 
the Qjot — Po overlay each other. Such an ambiguous situ- 
ation leads to the question which value of ai n j one should 
really assume for a given source to apply the presented 
method of the dynamical age estimate. This issue will be 
addressed further in section 4.3 below. 



4.2. Dependence on (3, AR, and 9 

The fitting results indicate that the effect due to changing 
particular values of the (3, AR and 6 parameters is of minor 
importance for the resulting age solution. Fig. 5 shows the 
influence of the exponent [3 in the density distribution of 
the ambient medium and of the source axial ratio AR = 
Rt x 2 on the age estimate. (Though the AR parameter is 
the observed one, it is much less precisely measured than 
the source linear size D.). The influence of the unknown 
inclination angle 9 on the age is also small. Comparison 
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Fig. 3. The influence of the initial energy distribution pa- 
rameters on the 'goodness' of fit quantified by the A mea- 
sure for CygA. 
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of Eqs.(l) and (3) yields t oc (D/ sin0)( 5 -W 3 . Thus, t 
increases no more than 1.21 and 1.15 for (3=1 and /3=2, 
respectively, if 9 decreases from 90° to 60°. We conclude 
that the influence of 9 on t is comparable to that caused 
by the uncertainty in the AR determination. 

4.3. Dependence of the age solution on kinetic energy 
delivered by the jet 

As noted in section 4.1, the fitted jet power increases with 
decreasing age (and hence with increasing «i n j ) for a given 
size and brightness of the source (see Figs. 1 and 2). In 
Fig. 6 we show both the fitted jet power and the age of 



8 



J. Machalski et al.: Dynamical age of FR II-type radio sources 



30 



CygA 



20 



10 



1 1 1 1 1 1 1 1 


iii 

7i,max=10 4 




1 A 10 - 
7i,max=10 




^inj 0.5 ( 

9 = 80° 


^^^^ 








0.6 1.0 1.4 1.0 

R T or , 



2.2 



2.6 



3.0 



Fig. 5. The influence of /3 and Rt on the age estimate for 
CygA, with a set of the fixed model parameters. 



42 



CygA 



Of 



M 
O 



38 
0.35 





\ 


7»ax-10 4 " 


\ 


\ i 


7max=10 


-16 V« 


sA Qjctt 




- 1 2 J7 






A-- 


\\ a*"/ 




- 6 





1 



54.4 



53.9 



53.4 



M 
O 



0.45 



0.55 



0.65 

a inj 



0.75 



0.85 



52.9 



52.4 
0.95 



Fig. 6. Diagram of the jet power Qjet, the age t, and the 
total energy Qj e t t delivered by the Cyg A jet to its co- 
coon, as functions of a\ n \ for the two (extreme) values of 



7i,max- The arrows indicate ai n j values corresponding to 
the minimum total cocoon energy. 

Cyg A as a function of a- m j for the two extreme values of 
7i,max considered in this section. As the parameters Qj et 
and t depend on ai n j in the opposite way, there is always 
a minimum value of the product Q^ ct x t. This product 
gives the kinetic energy delivered to the cocoon by the 
jet. The fitting procedure indicates thus that this total 
energy has a minimum for some specific value of a m j and 
fixed observational parameters. 

Basing on the above effect, we consider the age corre- 
sponding to the minimum kinetic energy as the 'best solu- 
tion' for the dynamical age of a given source (see section 



2.4). The 'best solution' for CygA gives i=10.4±1.6 Myr 
(a inj =0.574 ± 0.011) and t=5.7 ± 0.4 Myr (a inj =0.621 ± 
0.007) for 7i. max =10 4 and 7i !m ax=10 10 , respectively. Since 
one should expect in reality 7i imax in between of the two 
extreme values considered here for illustrative purposes, 
we safely conclude that the presented method indicates 
the age and the injection spectral index of CygA to be 
t = 6 — 10 Myr and a- m j = 0.57 — 0.62, respectively, in 
a very good agreement with other dynamical and spec- 
tral estimates presented in the literature (see, e.g., Kino 
& Kawakatu 2005). In addition, the corresponding kinetic 
power of the CygA jet is Q jct = (1 - 2) x 10 39 W, which 
is a factor of 2 — 3 higher than the value estimated by 
K2000, and a factor of 2 lower than the one obtained by 
Kino & Takahara (2004) by means of detailed analysis of 
the CygA hotspots' radiative properties. 

5. Fitting results for other radio galaxies 

5.1. The 'best solution' for the age 

Further results of the age estimation derived with our 
method are presented in more details for two other 
sources: giant radio galaxy B 1312+698 and normal-sized 
galaxy 3C55. The appropriate <5j ot vs. po diagrams for 
the given values of a; n j and 7i imax are shown in Figs. 7a 
and 8a. As previously, the dotted diagonal lines indicate 
an age related to the values of Qj ct and po fitted to ob- 
servations at the four different frequencies. The ideal, in 
this case, intersections of the four curves gives the precise 
'age solutions', as well as the precise estimates of the jet 
powers and the central core densities for both sources. 

The dependence of these 'age solutions' (marked as t 
for simplicity) on a; n j and 7i, m ax is presented in Figs. 7b 
and 8b, together with the Qjct — Q4nj and <5j ot t — a; n j 
curves. Clearly, the 'best solutions' exist within the theo- 
retically acceptable range of 0.5 < a; n j < 0.7. 

The 'best solutions' for the dynamical age of all the 
sources from Table 2 with the corresponding values of a; n j 
are collected in Table 3. In order to show how these solu- 
tions depend on the assumed value of the maximum elec- 
tron energies, we give the best solution for 7i. max = 10 7 
and 3 x 10 8 . As shown, the differences are minor. In par- 
ticular, column 6 of Table 3 gives the ratio of the standard 
deviation to the mean of the ages derived with these two 
different values 7i. max , showing that the resulting age es- 
timates differ from ~4% to ^30%, mostly by ~10% only. 



5.2. Comparison of the 'best solution' dynamical ages 
with the synchrotron ages 

A quality of the dynamical age estimated with the method 
presented here can be justified by its comparison with the 
synchrotron (radiative) age, r rac j, derived from the spec- 
tral ageing analysis (e.g. Myers & Spanglcr 1985; Carilli 
et al. 1991). In this analysis the radiative age of emitting 
particles is determined from the 'break' frequency in the 
observed radio spectrum and the magnetic field strength, 
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usually computed under the equipartition conditions. In 
order to find this break frequency, the a; n j value must 
be known. In most of the published papers concerning the 
above analysis, the value of a; n j was usually identified with 
the slope of the observed low-frequency spectrum. Because 
of a reasonable criticism of this simple method (see sec- 
tion 1), Murgia (1996) developed the software SYNAGE 
which allows the best fit of the spectral data to the theo- 
retical models of the energy losses by constraining the val- 
ues of the most important model parameters, especially a 
value of ami- 



The ainj value for the C.I. model, aci, fitted to the 
spectrum of the lobes of the sources from Table 2 using 
the SYNAGE algorithm, is given in column 7 of Table 3. 
It is easy to see that the values of aci are typically higher 
than the values of a; n j determined by minimizing the total 
energy as proposed in this paper, though the errors in «i n j 
are of the same order as the errors in aci (also given in 
column 7). Only for CygA and (formally) for 3C165 the 
fitted values of aci are lower than the values of ai„j . 

A discrepancy between the values of aci and ai n j for 
the other sources can be explained either by (i) a lack of 
flux density data at frequencies below ^100 MHz (which 
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Table 3. a in j values relating to the minimum total energy of the jet for two different values of 7i !max , and the 
corresponding 'best solution' of age for the analyzed sample of FR II-type radio galaxies. Column 6 gives the ratio of 
the standard deviation to the mean of the ages in columns 3 and 5. Columns 7 and 8 give a; n j value fitted for the C.I. 
model using the SYNAGE package of Murgia (1996), and the synchrotron (radiative) age of the sources compiled 
from the literature, respectively 
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may be crucial for a proper determination of the slope of 
the initial spectrum) and the uncertainties of the available 
data, or by (ii) an influence of other physical processes not 
taken into account in the KDA model, such as reacceler- 
ation of the relativistic particles, non-standard evolution 
of the cocoon's magnetic fields, etc. Note, that the aci 
values, fitted with the SYNAGE and being much higher 
than the effective values of a- m j determined by the min- 
imum jet energy delivered to the source (cocoon), may 
lead to a most unlikely young age and a corresponding 
enormously high jet-head advance velocity, iv For exam- 
ple, ainj ~ aci=0.898 for the high-rcdshift radio galaxy 
3C294 (see Table 3) would yield the dynamical age of 0.25 
Myr only, which would imply the unacceptable high value 
of i> h =0.936c! 

We note in this context, that large (~ 0.3 c) advance 
velocities are, in fact, observed, however only in 'compact 
symmetric objects' (CSOs; see, e.g., Owsianik & Conway 
1998, Polatidis & Conway 2003, Gugliucci et al. 2005). 
These sources, characterized by the projected linear sizes 
less than ~1 kpc, arc believed to be young versions of the 
extended radio galaxies, possibly even of the powerful clas- 
sical doubles (see a review by O'Dea 1998). And indeed, 
the implied kinetic ages (between 100 yr and 10 4 yr) usu- 
ally agree well with the spectral ages determined for some 
of these objects assuming energy equipartition within the 
lobes (Katz-Stone & Rudnick 1997). On the other hand, it 
should not be surprising that the advanced velocities of the 
jets in CSOs are larger than the ones in 'regular' FR II ra- 
dio galaxies, for which Scheuer (1995) gives an upper limit 
of 0.1 c. That is because of the differences in the lobes/jets 
evolution for compact and extended objects, caused for ex- 
ample by the differences in the ambient medium density 



profiles. In particular, Kawakatu & Kino (2006) argued for 
constant jet advanced velocity of evolving young/compact 
radio sources, while for the decreasing advanced velocity 
with increasing linear size at later evolutionary stages (see 
also in this context Alexander 2000). Therefore, here we 
take 0.3 c as an upper limit for the 'realistic' values of the 
jet advanced velocities in all the FR II radio galaxies con- 
sidered in this paper, in agreement with Scheuer (1995), 
and multi-epoch observations of young radio sources. 



Finally, the synchrotron ages r ra d (with associated er- 
rors) for all the sources in Table 2 — collected from the 
papers of Alexander & Leahy (1987), Leahy et al. (1989), 
Liu et al. (1992), and Parma et al. (1999) — are given in 
column 8 of Table 3. As shown, the synchrotron ages and 
our dynamical age estimates are comparable for the young 
sources, i.e. those whose ages are less or much less than 
10 Myr. Such behavior is in fact consistent with the con- 
clusions of Blundell & Rawlings (2000), who argued that 
the spectral ageing analysis is most likely erroneous for 
lobes older than 10 Myr. Indeed, in the case of the largest 
(and oldest) sources from Table 3 (3C 55, 3C 103, 3C 292, 
B1209+745, B1312+698), the ratios of the dynamical age 
estimates to the synchrotron ages (see column 9 of Table 
3) rise to a factor of 2 — 4. This agrees also with the ear- 
lier result on this ratio determined for the three giant ra- 
dio galaxies: J0912+3510, J1343+3758, and J1451+3357 
(Jamrozy & Machalski 2005). The 'best solution' for the 
age of the analyzed radio galaxies vs. their synchrotron 
age is plotted in Fig. 9. 
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Fig. 9. 'Best solution' age vs. the synchrotron age for the 
radio galaxies from Table 2. 

5.3. The jet power, the central core density, and other 
physical parameters 

Given the age determined for a given source, we can derive 
its other physical parameters. It was shown in section 4.1 
that each 'age solution' determined by the intersection of 
the age curves in the Qjet — Po diagrams (Figs. 2, 7a, and 
8a) indicates unique values of the jet power Qj et and the 
central core density po- It also provides unique values of 
the source (cocoon) energy density and pressure, u c and 
p c , respectively, as well as the total radiative energy E tot 
(see Table 1). We note, that the minimum total cocoon en- 
ergy (which gives the 'best solution' of the age) naturally 
corresponds to the equivalent minimum of u c and p c . 

All the physical parameters corresponding to the 'best 
solution' of age for the sources considered in this paper, 
obtained with 7i. max — 10 7 and given in Table 3, are col- 
lected in Table 4. The last column of Table 4 gives in 
addition the ratio of the mean expansion velocity of the 
jet head to the speed of light, i.e. one half of the projected 
linear size of a given source divided by its age estimate. 

5.4. A limitation of the method 

A justified question may be risen: starting from what 
age/size (or equivalently stage of the source's evolution) 
the described method can be applied? A priori the method 
is applicable to the sources for which radio emission is 
dominated by the lobes or cocoons evolving within a 
power-law ambient medium density profile. Therefore, our 



method cannot be applied directly to the young (CSO- 
type) sources whose linear sizes are < 1 kpc, because it 
is unlikely that their observed radio emission may be at- 
tributed entirely to the non-relativistic lobes. In particu- 
lar, in many of such sources contribution from the underly- 
ing relativistic jets and bright mildly-relativistic hotspots 
to the total radiative output may be significant. Finally, 
it is also not obvious if some (small) fraction of CSOs is 
'frustrated' by the dense ambient medium, as advocated 
by many authors (sec, e.g., Dc Young 1993, Alexander 
2000). 

Slightly larger versions of CSOs, called 'medium 
symmetric objects' (MSOs, or 'compact steep-spectrum' 
sources, CSS, if the spectral characteristics are consid- 
ered), are probably more promising in this context, since 
their projected linear sizes are typically 1 — 10 kpc. For ex- 
ample, the dynamical age estimates determined with our 
method for two of such objects from the sample of Murgia 
et al. (1999) are 142±35 kyr (a inj =0.571±0.020) for 3C 67 
with D = 11 kpc, and 258±67 kyr (a inj =0.753±0.100) for 
3C 186 with D = 15.7 kpc. The above dynamical age esti- 
mates are larger by a factor of 2 — 3 than the synchrotron 
ages found by Murgia et al. (which are 51 and 113 kyr 
for 3C 67 and 3C 186, respectively). However, such a dis- 
crepancy may be expected, as these synchrotron ages were 
determined using the break frequency in the spectrum of 
the entire source. That spectrum must be dominated by 
the radiation produced by the hottest regions of a given 
object, likely very close to a jet terminal shock. For this 
reason, the spectral ages found by Murgia et al. can be 
significantly lower than the dynamical ages derived by us. 
On the other hand, it is worth to emphasize that the ai n j 
values fitted by Murgia et al. in the frame of the C.I. model 
and our values determined from the minimum jet kinetic 
energy agree between themselves very well. 

6. Remarks on the injection electron spectrum 

Jets in powerful radio sources were considered as being 
characterized by sub- or eventually mildly-relativistic bulk 
velocities on large (hundreds of kpc) scales, Vj ct < 0.95 c 
(e.g., Wardle & Aarons 1997). Also, advanced velocities 
of their terminal features, the hotspots, were known to be 
non-relativistic, «h s < 0.1 c (Scheuer 1995). For these rea- 
sons, it was widely assumed that the double shock struc- 
tures formed at the heads of the jets in quasars and FR II 
radio galaxies were in fact always non-relativistic phenom- 
ena. In particular, it was argued that the non-relativistic 
reverse shock formed at the jet terminal point was the 
place of the efficient diffusive (Fermi I) particle accelera- 
tion process, leading to formation of the power-law par- 
ticle energy spectrum within the hotspot, injected later 
to the expanding lobes/cocoons by means of the plasma 
backflow. 

It is known that the well-established theory for parti- 
cle shock acceleration in the non-relativistic test-particle 
limit predicts a universal power-law index of the accel- 
erated particle continuum p = 2, weakly dependent on 
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Table 4. Physical parameters of the sources from Table 2 



Source 


log(Q je t[W]) 


logWkgm s ]) 


log(w c [Jm a ]) 


log(p c [Nm-^]) 


log(£ t ot[J]) 


Qjett/ Etot 


Vh/C 


CygA 


39.19 


-22.45 


-10.42 


-10.59 


53.17 


2.24 


0.032 


3C55 


39.28 


-22.02 


-11.24 


-11.42 


53.68 


3.38 


0.030 


3C103 


38.54 


-21.92 


-11.63 


-11.81 


53.10 


3.70 


0.017 


3C165 


37.99 


-22.46 


-12.06 


-12.22 


53.08 


2.06 


0.008 


3C239 


39.90 


-22.91 


-9.93 


-10.11 


53.60 


1.84 


0.054 


3C247 


39.19 


-23.75 


-10.50 


-10.68 


52.51 


2.61 


0.093 


3C280 


39.52 


-23.50 


-10.51 


-10.68 


53.17 


1.90 


0.071 


3C292 


39.18 


-22.29 


-11.83 


-12.01 


53.63 


4.54 


0.082 


3C294 


39.95 


-22.65 


-10.01 


-10.18 


53.67 


2.21 


0.060 


3C322 


39.91 


-23.30 


-10.82 


-11.00 


53.58 


2.70 


0.115 


3C330 


39.41 


-22.71 


-11.15 


-11.33 


52.93 


6.26 


0.098 


3C332 


37.45 


-22.61 


-12.06 


-12.24 


52.08 


2.67 


0.010 


B0908+376 


36.52 


-23.62 


-12.23 


-12.41 


51.06 


1.73 


0.006 


B1209+745 


36.98 


-22.80 


-13.37 


-13.54 


52.17 


3.07 


0.009 


B1312+698 


37.20 


-23.00 


-13.25 


-13.42 


52.33 


2.93 


0.012 



the plasma conditions on both sides of the shock (e.g., 
Blandford & Eichler 1987). Interestingly, radio spectral 
indices of powerful hotspots are indeed close to — though 
not always and not exactly — the corresponding value 
a inj = {p — l)/2 = 0.5. Such an agreement supported 
the 'non-relativistic shock' paradigm for the origin of non- 
thermal electrons within heads and cocoons of FR II type 
radio sources. It stimulated further development of this 
idea in the context of modeling the non-thermal hotspots' 
emission (Heavens & Meisenheimer 1987), enabling the 
successful explanation of several observational features 
(e.g., Meisenheimer et al. 1989, 1997). An important im- 
plication is that in a framework of such a 'standard' model 
the injection spectral index for all the cocoons of powerful 
radio sources is expected to be strictly a- m j — 0.5. 



As mentioned in section 2.3, a number of low-frequency 
radio observations of FR II sources indicate that there 
is no universal injection spectral index for their cocoons 
(see in this context Rudnick et al. 1994). This conclusion 
is supported also by the detailed spectral ageing studies 
of the few brightest objects of this type, such as CygA 
(Carilli et al. 1991). In addition, the standard assump- 
tion regarding non-relativistic velocities of the large-scale 
jets was recently put into question by the most recent X- 
ray observations (e.g., Celotti et al. 2001, Tavecchio et al. 
2005). If this is indeed the case, then one should expect 
significant modifications in the energy spectrum of the ac- 
celerated particles when compared to the non-relativistic 
limit. Moreover, as argued below, also the assumption re- 
garding a single power-law form of the injected electrons 
is most likely far from being realistic. For these reasons, 
in our analysis we decided to treat basically free 

'effective' parameter, which may be different for different 
sources, and which may be far from the canonical value of 
0.5. 



6.1. Low-energy electrons 

In this paper we assumed that the power-law continuum of 
electrons injected by the hotspot to the cocoons extends 
down to electron Lorentz factors 7i, m i n = 1. In general, 
synchrotron emission of the electrons with 7; < 10 3 can- 
not be observed directly from the hotspots, since it corre- 
sponds to frequencies < 100 MHz for a typical value of the 
hotspots magnetic field ~ 100 /iG. Thus, the assumption 
regarding the low-energy electron spectrum cannot be ver- 
ified by observations in a simple way, and should rather 
be constrained by the theory of particle acceleration. 

Unfortunately, such a theory is not yet well developed. 
We note in this context, that powerful jets at large scales 
are most likely dynamically dominated by cold protons 
(see, e.g., Sikora et al. 2005 and references therein). This 
means that the thickness of the terminal reverse shock 
should be roughly a few gyroradii of the non- or mildly- 
relativistic protons, and hence that the energies of the 
electrons which are able to undergo diffusive shock ac- 
celeration must already be ultrarelativistic, 7; > m p /m c . 
The questions are then (i) which process is responsible for 
acceleration of 1 < 7; < 10 3 — 10 4 electrons, (ii) is the 
resulting low-energy part of the electron continuum char- 
acterized by a power-law form at all, and if so (iii) is the 
low-energy spectral index the same as for the high-energy 
part of the electron continuum (71 > 10 3 — 10 4 ), originat- 
ing (by assumption) from the diffusive shock acceleration 
process. 

Such problems in the context of FR II hotspots were 
not discussed much in the literature, with only a few pa- 
pers addressing these important issues. For example, Kino 
& Takahara (2004) analyzed limiting possibilities for the 
efficient and non-efficient temperature coupling between 
dynamically dominating protons and electrons at the jet 
terminal shock, resulting in the formation of an ultrarel- 
ativistic maxwcllian-type electron distribution at energies 
7i ~ r jot m p /m e and 7; ~ r jet , respectively (where r jot 
is the jet bulk Lorentz factor), followed by the shock- 
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formed high-energy power-law tail. They argued, that 
multiwavelength observations of the CygA hotspots con- 
tradict both possibilities. This implies either an electron- 
positron jet composition (as preferred by Kino & Takahara 
2004) , or simply more complex (broad-band) energy spec- 
trum of low-energy electrons (instead of an almost mo- 
noenergetical maxwellian-like distribution). In fact, for- 
mation of such a broad-band energy spectrum seems to 
be more likely, since in the case of extragalactic jets any 
energy/temperature coupling between protons and elec- 
trons cannot be collisional in nature, but must be medi- 
ated by the magnetic field and plasma turbulence (see in 
this context Hoshino et al. 1992). 

To sum up, if the dynamical domination of the cold 
protons within powerful extragalactic jets is the case, one 
should expect (possibly significant) differences in the spec- 
tral shape of the low-energy (1 < 7; < 10 3 — 10 4 ) and high- 
energy (7; > 10 4 ) parts of the electron continuum formed 
at jet terminal shocks. Indeed, low-frequency flattening 
observed in the radio emission of CygA hotspots (Carilli 
et al. 1991), as well as the most recent analysis of the 
X-ray emission of the lobes of giant radio galaxy 4C 39.24 
(Blundcll et al. 2006) , are both consistent with such an ex- 
pectation. This in turn implies that the broad-band elec- 
tron spectrum injected by the hotspots to the expanding 
cocoons is in reality not of a simple single power-law form. 

6.2. Relativistic shock velocities 

In the section above we argued that the low-energy part 
of the electron distribution formed at the hotspots of 
powerful jets is most probably broad-band, possibly even 
power-law, but with the (unconstrained at the moment) 
power-law slope most likely different from the one charac- 
terizing its high-energy part. The difference may be due 
to the fact that while the low-energy part is formed by 
poorly known processes of energy exchange between differ- 
ent plasma species within the colissionless shock front, the 
high-energy part is supposed to originate from the 'stan- 
dard' Fermi I acceleration mechanism. Yet another issue 
is if this high-energy part always possesses some universal 
value of its power-law slope, the same for every powerful 
radio galaxy. We believe not, since even if the diffusive 
shock acceleration process is indeed responsible for for- 
mation of the latter, the expected relativistic effects will 
modify the resulting electron spectrum when compared to 
the non-relativistic test-particle model value p — 2. 

As already mentioned, there is growing evidence that 
relativistic velocities are present in large-scale FR II and 
quasar jets, and so the reverse shocks formed at the jet 
heads are at least mildly relativistic. An important fact in 
this context is that polarization studies usually indicate 
magnetic fields perpendicular to the jet axis within the 
hotspots (even in the case of giant-sized radio galaxies, see 
Machalski et al. 2006), and hence an oblique orientation of 
the mean magnetic field lines with respect to the reverse 
shock front. It was noted by Begelman & Kirk (1990), 



that oblique relativistic shocks are likely to be superlumi- 
nal, in a sense that the shock velocity projected on the 
magnetic field vector exceeds speed of light. In such a sit- 
uation, diffusive (Fermi I) acceleration process may not be 
possible at all, since particles moving on average along the 
magnetic field lines are not able to cross the shock front 
more than once. Moreover, as shown recently by Niemiec 
& Ostrowski (2004), even if the magnetic field configura- 
tion is luckily such that the diffusive particle scattering at 
both sides of the (mildly-)relativistic shock is in principle 
possible, the resulting energy spectrum of the accelerated 
particles strongly depends on the turbulence conditions 
in vicinity of the shock. As a result, a variety of spectral 
shapes for the accelerated particles is obtained, including 
different spectral indices and deviations from power-law 
forms. We believe that this may be partly the reason for 
different values of the injection spectral index obtained 
here for different sources. 

6.3. Jet/hotspot intermittency 

Multiple hotspots detected at the same sides of many 
FR II radio galaxies and quasars indicates that the jet 
terminal features are intermittent in nature (see, e.g., 
Scheuer 1982, Williams & Gull 1985, Leahy & Perley 
1995). Typically, a pair of the hotspots is observed at 
one side of the cocoon, with one feature being stronger 
and more extended ('secondary hotspot'), and the other 
one being weaker and more compact ('primary hotspot'), 
as in the famous source CygA. The secondary and pri- 
mary hotspots are identified with 'older (extinguishing)' 
and 'younger (recent)' jet termination points. We note, 
that such an intermittency may be a direct consequence of 
the intermittent nature of the jet phenomenon, discussed 
in the literature for a number of reasons (e.g., Reynolds 
& Begelman 1997, Stawarz et al. 2004). 

Whatever the case is, the issue of interest here is 
that multiple hotspots differ in multifrcquency radiative 
and geometrical properties (see Hardcastle et al. 1998). 
This should in fact be expected as a natural consequence 
of the ageing/evolutionary processes taking place in a 
mildly-relativistic shock front 'fed' by intermittent jet 
outflow. An important consequence of this fact is that 
the electron energy distribution injected to the cocoons 
by the hotspots, characterized by lifetimes much shorter 
than the lifetime of the whole radio source, is in fact an 
average over different activity phases of the jet terminal 
features. As a result, the injection spectral index inferred 
from the analysis of the cocoon's radio continuum may 
differ substantially from the low-frequency spectral index 
of the currently active hotspot in a given radio source, as 
suggested in our modeling. 

The above discussion, pointing out all the reasons 
why one should not expect a single broad-band power-law 
shape of the electrons accelerated at the jet head and 
injected in to the expanding cocoons, agrees with the 
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results of detailed investigation of radio continuum in 
FR II sources by, e.g., Rudnick et al. (1994), Treichel et 
al. (2001), and Young et al. (2005). These authors, by 
applying sophisticated technique of radio data analysis, 
found that there is no universal injection power-law 
spectral index for the cocoons of powerful radio sources, 
that the injected population of relativistic electrons 
thereby is characterized by complex ('curved') spectral 
shape, and that there exists a variety of distinct spectral 
features within the extended lobes. We argue, however, 
that for the purpose of simple estimates regarding the 
sources' main parameters (such as ages, jet powers, etc.) 
as presented here, a single power-law injection spectrum 
may be a relatively good approximation. The effective 
spectral index of such a model particle energy distribution 
may however be different for different sources. In this 
context we note that somewhat surprisingly our analysis 
suggests a relatively narrow range for this parameter, 
namely p = 2.0 - 2.4. 

7. Conclusions 

Here we propose a new approach to determining the ages 
of powerful radio sources, on the one hand exploiting a 
dynamical model developed for these objects by Kaiser 
et al. (1997; KDA), and on the other hand using multi- 
frequency radio observations not necessarily restricted to 
the high-resolution ones. The KDA model is chosen for its 
minimal number of free parameters when compared to the 
more sophisticated scenarios by Blundell et al. (1999) or 
Manolakou & Kirk (2002). 

We apply the assumed dynamical model to a number 
of FR II type radio galaxies observed at different radio 
frequencies, and fit — for each frequency separately - 
the model free parameters to the sources' observed quan- 
tities. Such a procedure, enlarging in fact a number of 
observables, enables us to determine, relatively precisely, 
ages and other crucial characteristics (such as the jet ki- 
netic power) of the analyzed sources. We show, that the 
resulting age estimates extend from about 1 Myr for the 
small-size but powerful high-redshift radio galaxies to over 
100 Myr for the giant-sized radio galaxies. The estimated 
ages agree very well with the ones obtained by means of 
'classical' spectral ageing method for objects not older 
than 10 Myr, for which good-quality data are available. 
Meanwhile, in the case of sources older than 10 Myr, our 
age estimates are systematically larger (by a factor of 2—4) 
than the ages obtained by means of spectral ageing anal- 
ysis. Such a behavior is in fact consistent with the conclu- 
sions of Blundell & Rawlings (2000), who argued that the 
spectral ageing analysis is most likely erroneous for the 
lobes older than 10 Myr. 

Interestingly, the ages of FRII-type sources estimated 
by using the presented method correspond always to 're- 
alistic' values of the jet advance velocity < 0.3 c, i.e. con- 
sistent with those measured directly with the VLBI tech- 
nique in many of the CSO-type sources (see section 5.2.). 
Thus, we believe that the method presented here is appli- 



cable also in the case of the oldest radio galaxies, and/or 
ones for which the only available low-resolution radio data 
do not allow for detailed spectral ageing studies (i.e., for 
majority of powerful radio sources). This implies that the 
applied procedure of favoring the value of the injection 
spectral index which corresponds to the minimum of the 
cocoons' total energy (see section 4.3), although to some 
level arbitrary, nevertheless gives very reasonable results. 

Our analysis indicates that the main factor precluding 
precise age determination for FR II type radio galaxies 
regards the poorly known shape of the initial electron en- 
ergy distribution injected by the jet terminal shocks to the 
expanding lobes/cocoons. We consider briefly this issue, 
and conclude that the broad-band single power-law form 
assumed here (and in all the other analogous models) may 
be accurate enough for the presented estimates, although 
most likely it does not strictly correspond to some well- 
defined realistic particle acceleration process. Instead, it 
should be considered as a simplest model approximation 
of the initial electron continuum, averaged over a very 
broad energy range and over the age of the source. The 
effective spectral index of such a distribution may be dif- 
ferent for different sources, however within the relatively 
narrow range p = 2.0 — 2.4 suggested by our modeling. 

In a subsequent paper we indend to apply the pre- 
sented method for the age determination in a large sample 
of FR II-type (but also younger, MSO-type) radio sources. 
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